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Abstract 

The behaviour of a driven double well Duffing-van der Pol (DVP) oscillator 
for a specific parametric choice (| a |= /?) is studied. The existence of dif- 
ferent attractors in the system parameters (/ — uj) domain is examined and 
a detailed account of various steady states for fixed damping is presented. 
Transition from quasiperiodic to periodic motion through chaotic oscillations 
is reported. The intervening chaotic regime is further shown to possess is- 
lands of phase-locked states and periodic windows (including period doubling 
regions), boundary crisis, all the three classes of intermittencies, and transient 
chaos. We also observe the existence of local-global bifurcation of intermittent 
catastrophe type and global bifurcation of blue-sky catastrophe type during 
transition from quasiperiodic to periodic solutions. Using a perturbative pe- 
riodic solution, an investigation of the various forms of instablities allows one 
to predict Neimark instablity in the (/ — iv) plane and eventually results in 
the approximate predictive criteria for the chaotic region. 
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I. INTRODUCTION 



The Duffing-van der Pol oscillator (DVP) 

X — //(I — x^)x + ax + (3x^ = f cosut, > 0, ^— = (1) 

is a ubiquitous nonlinear differential equation which makes its presence in physical, engi- 
neering and even biological problems [1-7] . It is a generalization of the classic van der Pol 
oscillator equation. It can be considered in atleast three physically interesting situations, 
wherein the potential V{x) = ^ + ^ is a (i) single well {a > 0,P > 0), (ii) double well 
(a < 0, /5 > 0) or a (iii) double hump (a > 0,/3 < 0). Each one of the above three cases 
has become a classic central model to describe inherently nonlinear phenomena, exhibiting 
a rich and baffling variety of regular and chaotic motions. 

Chaotic motion in system (1) with a single well type restoring force was investigated by 
Ueda and Akamatsu [8] as a model of negative resistance oscillator and later on was studied 
by a number of other authors [9-12], who noted symmetry breaking of attractors and the 
onset of chaotic dynamics. Bountis etal.[ll] have investigated the nonintegrability of a family 
of DVP oscillators by studying analyticity properties of the solution in the complex time 
plane and proved that infinitely sheeted structure (ISS) exists in this system. Rajasekar, 
Parthasarthy and Lakshmanan [13] pointed out that the DVP oscillator with double well 
potential exhibits Smale-horse shoe chaos when transverse intersections of the homoclinic 
orbits occur. Further Kao and Wang [14] had analog simulated the DVP oscillator with a 
double hump potential and discussed the various mode locking, multiple hysteresis, period 
doubling route to chaos, intermittent hopping and crises phenomena. 

Recently Szemplinska-Stupnika and Rudowski [12] reported that a single well type DVP 
oscillator exhibits chaotic motion between two types of regular motion, namely periodic 
and quasi periodic oscillations, in the prinicipal resonance region for a specific values of 
the parameter /(= 1.0) and a range of u values (0.8 - 1.0). Also they have obtained a 
perturbative solution for the periodic oscillation and carried out stability analysis of such 
solution to predict Neimark bifurcation. However no such analysis exists for the important 
case of double well type DVP oscillator so far in the literature, which is atypical in the sense 
that even in the absence of forcing it shows the existence of multiple attractors [6,7]. 

Considering the DVP oscillator with a double well type restoring force in the form 

X — /i(l — x'^)x— I q; I a; -|- (3x^ = f cos cot, (3 > (2) 

we notice that the three equilibrium points of the system (2) for / = correspond to 
— I q; I a; -I- f3x^ — 0, so that we have the stable fixed points a: = ^V^'^ ^^<i unstable 
fixed point which is hyperbolic at x^"^ = 0. Actually Xi 2 are elliptic points for \ a \— P 
and become stable foci for | a |> while they are unstable foci for | a |< As a result, 
the system (2) exhibits a large orbit (LO) motion, which always encircles all the three 
equilibrium points for the case \ a \ = (3. As far as | a |> /? is concerned, the system exhibits 
both small orbit (SO), that is oscillation around any one of the stable fixed points, and 
large orbit motion (LO), depending upon the values of the other control parameters and 
also initial conditions . 
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In this paper we undertake an investigation of the dynamics of the double-well DVP 
oscillator (2) and show that it is a rich dynamical system, possessing a vast number of 
regular and chaotic steady states. In particular, considering the special case \ a \= [3 (the 
case \ (3 is even more richer and the results will be presented separately [16]), we bring 
out the existence of transition from quasiperiodic to periodic motion in the (f—co) parameter 
space via chaotic motion . The novel features we observe are that in the chaotic sea there are 
many isles of periodic and phase-locked states, which exhibit period doubling phenomenon, 
intermittencies, crises etc., along with regions of transient chaos, corresponding to local 
bifurcations. Besides, there are also transitions from quasiperiodicity to period T orbits 
corresponding to global bifurcation of blue-sky and local-global bifurcation of intermittent 
catastrophes. We also present a perturbative approach to the study of bifurcations which 
occur near the principal resonance. The analysis allows us to derive the algebraic equations 
for the instability boundaries. 

The plan of the paper is as follows. In sec. II we present the numerical results for 
different steady states, bifurcation routes and chaos for system (2). In sec. Ill we develop 
a perturbative solution and obtain expressions for the stability regions and compare them 
with numerical results. Then, we compare the results with the dynamics of the double well 
Duffing oscillator in sec. IV. Finally, sec. V summarizes our results. 

II. NUMERICAL RESULTS 

Equation (2) is numerically integrated using the fourth order Runge-Kutta algorithm 
with adaptive step size with parameter values fixed at | a |=0.5, P —0.5, and /i—O.l inorder 

to study the large orbit behaviour mentioned in the introduction. The transitions are also 
characterized by tracing the time evolutions, phase portrait, Poincare map, Fourier spectrum 
and Lyapunov exponents. For identifing different steady states, the dynamical transitions 
are traced out by two different scanning procedures: (1) varying a; at a fixed / (frequency 
scanning) and (2) varying / at a fixed uj (amplitude scanning). The resulting phase diagram 
in the (/ — oj) parameter is shown in Fig.l. The diagram covers the transition thresholds 
in the region of principal and super harmonic resonances, 0.4 < c<j < 1.0 and the forcing 
strength lying in the region 0.0 < / < 0.2. The various features in the phase diagram are 
summarised and the dynamical transitions of the attractors are elucidated in the following. 

A. Phenomena of steady states 

One observes that equation (2) admits the free-running solution when the external force 
is absent (/ = 0.0). When it is present and for low / values and low uj values, the frequency 
of the system becomes incommensurate with the external frequency. Consequently, the 
system exhibits multifrequency quasiperiodic oscillations. When the value of the external 
frequency to exceeds certain critical value for fixed /, a transition from quasiperiodic to 
periodic oscillations occur on increasing u (see Fig.l) essentially due to supercritical Neimark 
bifurcation (see Sec. II B below and Sec III). This phenomenon continues until a critical / 
value (/ - 0.115). 
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On a further increase of the forcing parameter f,f> 0.115, the system exhibits chaotic 
motion between the two regular motions, that is quasiperiodic and periodic oscillations, 
within a range of the driving frequency u. For example, at / = 0.13, chaotic motion occurs 
in the region cu e (0.546,0.553) (see Fig.2a). 

As the forcing parameter / increases further, within a very narrow frequency region, 
chaos-periodic windows-chaos type transition is found to occur between the two regular 
oscillations. For example at /=0.14, period 5T solution occurs in the frequency range 
cu e (0.545,0.551) within the chaotic range lu G (0.525,0.553) (sec Fig.2b). At /=0.17, 
period doubling phenomenon occurs in the window region of the frequency u G (0.503, 0.529). 
In the above range, we note that the period 2 orbit is born at a; = 0.529;it undergoes a period 
2x2"^ doubling bifurcations, finally leading to the onset of chaos as u> decreases (see Fig. 2c). 

Higher values of the forcing strength / introduces the appearance of a new transition. 
When this increase is coupled with increasing frequency, the quasiperiodic motion suddenly 
changes into a phase locked attractor. As an example at / = 0.19, the transition from 
quasiperiodic oscillation to phase locked states of period 3 orbit occur at a; = 0.453 and 
this locked state persists in the frequency range uj G (0.453, 0.47), which is then followed by 
chaos, reverse period doubling, chaos and periodic solution (see Fig. 2d). 

The details of our numerical study are summarised in the (/ — to) phase diagrams, Fig.l 
and Fig. 2. It depicts the system parameter region where quasiperiodic, large periodic 
and chaotic attractors exist. The curves denoted QP(1),QP(2) and P(l) are the boundaries 
of transition from quasiperiodic to the chaotic state, quasiperiodic to periodic state and 
chaotic to periodic state respectively. Also one observes in the entire transition regions 
where coexistence of multiple attractors occur. Further, beyond the curve P(l), there are 
some regions which exhibit phase locked periodic and transient chaotic states. However, in 
this paper, those states are not discussed in detail. 

The various steady states as 1,2, 3,. ..,15 denoted in Fig. 1 are then illustrated in Fig.3. 
Regular attractors are illustrated by their phase portraits and quasiperiodic and chaotic 
attractors by their Poincare maps. 

The first three points (l)-(3) in Fig.3 are examples of almost periodic (quasiperiodic) 
orbits for low values of /. Then the points (4)-(15) are essentially located at the principal 
and super harmonic resonance regions for large values of the forcing parameter, / > 0.12, 
at increasing driving frequency u. We observe here the following: quasiperiodic orbit (point 
(4)), period 3T orbit (point (5)), chaotic orbits (points (6), (7)), period doubled orbits (points 
(8), (9)), chaotic orbits (points (10), (11)), period 5T orbit (point (12)) and period T orbits 
(points (13)-(15)). 

B. Classification of Bifurcations 

The complicated dynamical behaviours of the DVP oscillator (2) with | a |= /3 due to 
the presence of the double well restoring force has been confirmed by the phase diagram as 
discussed above. From the bifurcation theory point of view, these correspond to several types 
of bifurcations: secondary Hopf, intermittent and blue-sky catastrophes besides standard 
period doubling bifurcations which are discussed in the following section. 
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1. QP(1) region - Local bifurcations: Secondary Hopf (Neimark) bifurcation 



In analogy with the Hopf bifurcation, a bifurcation is expected at a critical value as the 
limit cycle loses its stability, so that an attracting torus is born. This is the secondary Hopf 
bifurcation or a Neimark bifurcation [17]. Further, the bifurcated solution can be either 
stable and supercritical or unstable and subcritical. For the present DVP oscillator (2) with 
I a 1= there is a very large transition region QP(1) corresponding to this secondary Hopf 
bifurcation. As an example, let us examine the transient process near / = 0.1. Figures 4 
show the Poincare maps with values of uj decreasing and with the starting values of x and 
X indicated as in brackets. In Fig.4a, for to — 0.59, we see a node like convergence to a 
point, while in Fig. 4b for u = 0.58 the convergence has a spiralling character and the rate of 
convergence is noticebly slower. In the last diagram of Fig.4c, at uj = 0.57, we see that the 
system is moving outwards from near unstable fixed point towards the attracting invariant 
closed curve and the approach is termed as supercritical. This is the typical behaviour for 
the curve QP(1) in Fig.l for / < 0.10. 

The other nature of Neimark bifurcation, namely subcritical behaviour has also been 
observed in system (2). For example, at / = 0.15, chaotic region exists in the range uj G 
(0.512,0.551) (see Fig.5). It follows that this narrow strip of chaotic motion is related to 
the transition from quasiperiodic to periodic oscillation via chaotic motion. This transition 
does not occur in a smooth, continous way as it is the case when the Neimark bifurcation 
is a supercritical one but occurs through a chaotic region. This type of transition has a 
close resemblence with the Duffing oscillator [15,17,18] case where it was shown that a lower 
frequency band of chaotic region is related to the saddle-node bifurcation, which causes 
a sudden change from/to T-periodic orbit to/from chaotic attractor and transient motions 
separate the two different steady states. Thus by analogy the occurrence of chaotic motion in 
the region of driving frequency that separates quasiperiodic and periodic oscillations can be 
interpreted as a subcritical Neimark bifurcation. This type of bifurcation has been reported 
earlier in ref[12] for the single well DVP oscillator, but the unstable motion corresponds to 
randomly transitional motion, whereas in the present double- well case this corresponds to a 
fully chaotic attractor and periodic windows. 

2. QP(2) region - Local-global bifurcation: Intermittent catastrophe 

Figs. 6 show the typical Poincare maps with forcing values / = 0.00172 and 0.00173 at 
UJ = 0.83 in the QP(2) region, where a transition from quasiperiodic to periodic motion 
occurs. In order to lock the quasiperiodic motion to the period T motion, the Poincare map 
points form a closed loop with points progressing more rapidly near the top of the attractor 
and more slowly where points are visibly dense. As the external force strength is increased, 
a pair of saddle and node develops such that quasiperiodic motion rapidly shrinks to the 
node, which is near the original attractor. This type of process is prototypical of intermittent 
catastrophe [17,19]. 
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3. QP(2) region - Global bifurcations: Blue-sky catastrophe 



As the frequency value is increased to a; = 0.998, a different transition process is captured 
for / = 0.05. Tfie typical Poincare map is shown in Figs. 7 with the values of uj at 0.998 and 
1. and / = 0.05. Prom Pigs. 7, we find that a periodic attractor bifurcates to a quasiperiodic 
attractor which is located inside the other, and the two attractors are typically disjoint 
and separated in phase space by a finite distance. Furthermore, it is not generic for a 
quasiperiodic attractor to appear suddenly at the same control threshold where periodic 
motion vanishes. The quasiperiodic attractor will have existed previously, or it will not 
exist at all; in either case the bifurcation will consist of periodic attractor simply losing 
stability - it vanishes into the blue. Such phenomenon is termed as blue-sky catastrophe 
[17,20] in the literature and this event involves collisions with saddle type objects. 

4- V-shaped region - Transitions to chaotic attractors 

We now enumerate the various possible attractors present in the system in the V-shaped 
region in Pig. 1. 

(i) Transient chaos and boundary crises: Chaotic behaviour is observed between 
boundary of curves QP(1) and P(l) in Pig.l. Near the boundary of each of the curves P(l), 
QP(1), the system behaves in a random way, with the trajectory moving in phase space as 
if it were on a strange attractor. However after a transient time, the motion settles into a 
regular attractor, that is near P(l) it settles into a periodic motion ( for example see Pig. 8) 
while near QP(1) into a quasiperiodic oscillation. Such a phenomenon is termed as transient 
chaos, which is a precursor to steady state chaos. Between these two transitional regions, 
periodic windows, phase locked states, chaotic attractor and period doubling phenomenon 
occur. In addition, the boundary crisis [21] of chaotic attractor appears as the value of the 
external frequency increases so that the dynamics corresponds to the curve P(l), where the 
boundary of the chaotic attractor touches the unstable periodic orbit. 

In addition, one observes the interesting fact that in the V-shaped region there exists 
different parametric values for which intermittency of all the three classic types occur. This 
seems to be rare in such low dimensional systems. 

(ii) Type I Intermittency: The parameter regions separating the periodic windows 
inside the V- shaped region correspond to various complicated dynamics including chaos. 
The precise stability boundary of each window has been found to correspond to a saddle-node 
instability. As an example, for / =0.17, if u! is increased across the saddle-node boundary, 
type I intermittency occurs[22-25]. One such intermittency signature is shown in Pig.9. The 
average laminar length (< Z >) of this type of intermittency is found to comply with the 
law < I >~ jjT^ with 5 ~ 0.52 ± 0.001 where ji — uj — u)'^ and w'^is the bifurcation threshold. 

(iii) Type II Intermittency: In the earlier section, we showed that two possibilities 
exist in the QP(1) region when the periodic motion encounters a Hopf-bifurcation. Either 
quasiperiodic motion results if the bifurcation is supercritical or a complicated evolution 
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appears if the bifurcation is subcritical. In the later case the transition to chaos is found 
to have intcrmittency signature in certain parametric regions. Such an intermittency sig- 
nature is shown in Fig. 10. A close look into the signature reveals that there are distinct 
phases of the regular motion which are punctuated by other phases which are apparently 
chaotic. According to the classical Pomeau-Manneville categorization of different types of 
intermittcncics based on local bifurcations, the present intermittency is of type II since the 
preceding bifurcation is a Hopf bifurcation [22-25]. The average length < / > of the laminar 
phase of this intermittency is found to comply with the law < / >~ {j^Y with 5=0.9321 
where /j, — cu — cu^ and co'^ is the bifurcation threshold. 

(iv) Type III Intermittency: Next we discuss yet another type of route inwhich the 

periodic orbits in the periodic windows are seen to undergo intermittent transition to chaos. 
One such intermittcny motion is shown in Fig. 11, which is a Poincare time series plot for 
/=0.14 and lj=0.53802. It is seen that the motion just before the onset of intermittency 
is of period 20 orbit which itself occurred due to the period doubhng of period 10 orbit. 
The laminar phase of Fig. 11 corresponds to period 40 orbit along with chaotic bursts. 
Therefore this is identified to have arisen out of a subcritical half subharmonic instability, 
that is subcritical period doubling. Thus according to PM classification this is type III 
intermittency [22]. The average length < / > of the laminar phase of this intermittency 
complies with the following scahng law predicted by Pomeau-Manneville: < I >~ {j^Y with 
S — 0.9912 where ji — uj — uj'^ and uj'^ is the bifurcation threshold. 



III. PERTURBATIVE ANALYSIS 

From the numerical studies reported in the earlier sections, we observed that the T 
periodic orbit within a range of driving frequency uj for low values of / is close to a harmonic 
function of time. Then by obtaining a first order approximate period T (^7^) solution and 
analysing the stability one can estimate the system parameters domain inwhich Neimark 
instability arises as was done in the case of single well DVP oscillator for fixed / and uj in 
ref [12]. Keeping this aim in mind wc look for a periodic solution of (2) using a perturbative 
method (with both | a \ and (3 fixed at 0.5). Applying the method of multiple scales [1,26] to 
equation (2), one can obtain the approximate solution about the stable fixed point Xg — 
in the form 

3 

X = + cicos{u!t + 4>) + — cos2{Lut + (p) iJ,sm2{ujt + 0), (3) 



where 



« = ^ , (4) 



^J{n^(a)-UJ''f + [l|Iu;a^] 



and fi^(a) = 1 — is the natural frequency of the autonomous conservative system (2) at 
/J, — and f — 0. 
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A. Linear stability analysis 



1. Soft-mode instability 

In order to examine the stability of the solution (3), we may look at a specific form 
of instability which manifests itself by an exponential growth with time of the harmonic 
components in the solution (3). This can be done by adding small disturbances to the 
amplitude and phase of the solution (3) as 

X ^ -^(a + SaY + (a + Sa) cos{u;t + (j) + 5(f))+ +^'^^) cos2((x;t + + 50) - 

^^^-^/xsin2(a;t + + 50). (6) 



Working out the linearized equation for 6a and 50, one ultimatively arrives at the following 
expression which corresponds to the first order instability limit as 



+ (^/x^a^ - 2a^ - 2) a;2 + (20^ + + l) = 0. (7) 



The above analysis is valid only for a fluctuation having the same frequency as the approx- 
imate solution x(t) considered. 



2. Hard-mode instability 



Now we examine another type of instability in which the perturbation may have dif- 
ferent harmonic components other than those in x(t). Following the spirit of the work of 
Szemplinska-Stupnicka and Rudowski [12], let us study the effect of a small disturbance to 
x{t) — x{t), where x{t) is the solution (3), in the form 



x{t) = x{t) + Sx{t). 
The linear variational equation for 5x{t) is then 

5x + Pi{t)5x + P2{t)6x = 0, 

where 



oo^ + y + ^ + TT 1 + 2aoacos6l-F ( ^ ) cos26'-F ... 



2 



a 
~2 



(8) 



(9) 



(10) 



P^(t) = - I a I -F3/3 ( ao^ + y + 2aoacos6'-F y COS26' 



+2/X 



q3\ ^3 

-ao^a a; sin ^ -|- aoO^to sin 29 — —uj sin 3d 



+ ... 



(11) 



and ao — — f a^; 02 = ^; 03 = 9 = out + (j), where dots correspond to higher harmonic 
components. Introducing now the transformation 
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Sx = u exp 



f 1 
— J P,{t)dt 



(12) 



equation (9) can be converted into a Hill's equation 

u + P{t)u^O, 

where 

p2 p 



(13) 



(14) 



Using the form of Pi and P2 given in (10) and (11), the transformation (12) can be rewritten 



as 



6x = u exp 



-At + 



COS 6* + y COS 29 j de 



(15) 



(16) 



Applying the Floquet theorem, one can look for a particular solution of (13) in the form 

M = exp(eii).0(i), (17) 

where (f){t) is a periodic function of time and ei is either real or imaginary. Thus the equation 
(15) becomes 



6x{t) = exp(ei - A)t.(j){t), 



where 



4){t) — (f){t) exp 



(200 a 



a 



cos9+— cos 26 1 d9 



(18) 



(19) 



The stability of solution (3) depends exclusively on the exponent coefficient (ei — A) in 
(18). Let now discuss the various possibilities to have stable solution. 

Case(i): Considering the case ei = ±iei so that ei is real and positive then 



Sx{t) = exp(±iei - A)t.4>{t). 



(20) 



It can be concluded that when ei is imaginary, the solution (3) is stable if A > or a > 
0.91 and unstable if A < 0. This form of instability (termed as Neimark instability) leads 
to a buildup of new harmonic components whose frequencies are incommensurate with the 

frequency of the periodic sohition (3). 

Case(ii): Considering the case ei = ±ei so that ei is real then 



5x{t) = exp(±ei - A)t.4){t). 



(21) 
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When ei is real, the solution (3) is stable if A > and > and this form of instability 
is approximately equal to the classic first order instability as given by Eq.(7). 

Therefore the form of instability defined by the equation (20) and so the condition 
A > leads to the build up of new harmonic components with frequencies uj + Ei and 
cu — ei. However, these frequencies are in general incommensurate with the frequency cu of 
the periodic solution (3) whereas for A < the solution is Tinstablc and so A = is the 
boundary of the instability. Thus, this instability can be interpreted as a Neimark instability, 
giving rise to a Neimark bifurcation. 

Let us now look at the resonance curves, shown in Fig. 12 and the two unstable regions 
defined by the condition (7) (first order instabihty) and the condition A < for Neimark 
instability. From the Fig. 12, the Neimark instability is expected to occur at the frequency 
value where the resonance curve crosses the critical boundary value a ~ 0.91. To determine 
the Neimark stability limit in the (/ — a;) parameter plane, we calculate the forcing pa- 
rameter / by using the resonance equation (4). Fig. 13 depicts the Neimark instabihty limit 
defined by the condition A < and the first order stabihty hmit described by the condi- 
tion (7). The numerical study results presented already in Fig.l are shown for comparison. 
The theoretically predicted Neimark instability values are reasonably close to the numerical 
results. 

IV. COMPARISON IN THE DYNAMICS WITH THE DOUBLE WELL DUFFING 

OSCILLATOR 

Finally it is of importance to compare the dynamics of the double- well Duffing oscillator 
as given by Szemplinska-Stupnika and Rudowski [15], 

X + jjLX— \ a \ X + [3x^ = f cos ut, (22) 

with that of double well DVP (2) for the same parametric values \ a \= j3 = 0.5 and 
H = 0.1. The results are compared in Table 1. In the lower frequency boundary region 
(see Fig.l of ref[15]), the Duffing oscillator exhibits symmetry breaking bifurcations while 
in the DVP, Neimark bifurcation boundary(Fig.l of the present paper) has been found. In 
the higher frequency boundary region, particularly in the prinicipal resonance region, the 
Duffing system exhibits period doubling bifurcations of small periodic orbits and cross well 
chaos. However in the DVP system, this region is always found to have highly regular orbits. 

Further, the Duffing oscillator exhibits period doubling route to chaos of large period T 
orbit in the lower frequency region, that is a; < 0.4, while in the present case, the system 
always exhibits almost-periodic oscillations in the region. In the present case, in every tran- 
sition boundaries such as QP(1), P(l) and QP(2), there are some regions where coexistence 
of multiple attractors arc found to occur. But in the Duffing oscillator case, the coexistence 
of multiple attractors are observed in the transition region from large period T orbit to 
cross well chaos. Blue-sky catasrophe, type II and HI intermittencies and various phase 
locked states are found to occur in the present case. However, in the Duffing oscillator case 
such phenomena have not been found (atleast to our knowledge). Naturally the dynamics 
exhibited by the DVP equation (2) is also quite distinct compared to the forced van der Pol 
oscillator [27]. 
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V. CONCLUSIONS 



Numerical studies show that the double well Duffing-van der Pol oscillator (2) with the 
parameter choice \ a\ = (3 exhibits a rich variety of attractors of periodic, quasiperiodic and 
chaotic types. Four varieties of transitions from quasiperiodic to periodic motions occur: (1) 
QP - periodic orbits (2) QP - chaos - periodic, (3) QP - chaos - periodic windows - chaos - 
periodic and (4) QP - phase locked states - chaos - periodic orbits. Besides these, local stable 
or supercritical, unstable or subcritical Neimark bifurcations and mode locking, intermittent 
catastrophe and blue-sky catastrophe bifurcations are also shown to exist. Transient chaos, 
period doubling phenomena, boundary crisis and intermittencies of all the three classic types 
are shown to occur and these were demonstrated with suitable examples in the (/ — uj) 
parameter space. In the literature, so far all the three intermittencies have been found to 
occur mostly in the higher dimensional or coupled systems [21-23]. However in the present 
case even in a single model of low-dimensional system, we are able to demonstrate all the 
three kinds of PM intermittencies. The various forms of instabilities of approximate periodic 
solution allows one to predict the Neimark bifurcation in the (/ — a;) parameter domain. 
Although some discrepency between true and theoretical predictions occurs, the approximate 
analysis throws some light to distinguish between the regular and chaotic regions. 
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FIGURES 



FIG. 1. Regions of different steady states exliibited by the double well DVP oscillator (2) at 
/i = 0.1. 

FIG. 2. (i) Bifurcation diagram for maximum amplitude x vs external forcing frequency a;, (ii) 
Maximal Lypanov exponet Xmax vs external forcing frequency uj of system (2). (a) QP - chaos - 
periodic orbit transitions for uj G (0.4,0.65) at /=0.13. (b) QP - chaos - periodic windows - chaos 
transitions for u G (0.4,0.65) at /=0.14. (c) QP - chaos - period doubling window - chaos transitions 
for u G (0.4,0.65) at /=0.17. (d) QP - phase locked states - chaos - period doubling windows - 
chaos transitions for u G (0.4,0.65) at /=0.19. 

FIG. 3. Various types of steady states at ^ = 0.1: (1)/=0.001, co =0.78; (2) /=0.001, u;=0.9; 
(3) /=0.12, uj=0A5; (4) /=0.15, uj=0A5; (5) /=0.19, uj=0A7; (6) 6j=0.515; (7) uj=0.52; (8) 
cc;=0.525; (9) a;=0.53; (10) /=0.15, w=0.54; (11) /=0.14, u;=0.53; (12) /=0.135, Cc;=0.54; (13) 
/=0.14, a;=0.58; (14) /=0.15, Cc;=0.58; (15) /=0.18, a;=0.58. 

FIG. 4. Trajectories near the Neimark bifurcation for system (2): /=0.1: (a) (j=0. 59(1,1); (b) 
a;=0.59(l,l); (c) a;=0.57(l,l). 

FIG. 5. (i) Bifurcation diagram for maximum amplitude x vs external forcing frequency uj, (ii) 
Maximal Lypanov cxponct Xmax vs external forcing frequency lo of system (2). QP-chaos-periodic 
orbit transitions for a; G (0.4,0.65) at /=0.15. 

FIG. 6. Poincare map of the system (2) before and after mode locking for (a) /=0. 00172, (b) 
/=0.00173, (c) /=0.00174 at Cc;=0.83. 

FIG. 7. Poincare map of blue-sky disappearence of Periodic attractor in system (2): (a) 
a;=0.998, (b) io=l. at /=0.05 



FIG. 8. Trajectories to show transient chaos for / = 0.125, a; = 0.535. 

FIG. 9. Signature of type I intermittency: Time series plot for /=0.17, (a) a;=0. 52601; (b) 
a;=0.526010001. 

FIG. 10. Signature of type II intermittency: Time series plot for /=0.12, (a) a;=0.5532; (b) 
a;=0.5537. 



FIG. 11. Signature of type III intermittency: Poincare time series plot for /=0.14, a;=0. 53802. 
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FIG. 12. Resonance curves and unstable regions of solution (3): I: branches unstable in the 
sense of first order instability; II : branches unstable in the sense of Neimark instability. 

FIG. 13. Regions of different steady states : Numerical (solid line) and theoretical (dashed 
line) stability limits. 
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TABLES 



I. Comparison of the orbits of Duffing and DVP oscillators 



Parameters 




Duffing Oscillator 


DVP Oscillator 




/ small 


small orbit 


QP orbit 


u g(0.4,0.6) 










/ large 


large orbit and chaos 


large orbit and chaos 


u G(0.6,1.0) 


/ small 
/ large 


small orbit 

small orbit and chaos 


QP orbit 
large orbit 
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